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Abstract. - We describe single electron tunneling through molecular structures under the 
influence of nano-mechanical excitations. We develop a full quantum mechanical model, which 
includes charging effects and dissipation, and apply it to the vibrating C60 single electron 
transistor experiment by Park et al. [Nature 407, 57 (2000)]. We find good agreement and 
argue vibrations to be essential to molecular electronic systems. We propose a mechanism to 
realize negative differential conductance using local bosonic excitations. 



Introduction. Experiments on electronic transport through nano-scale systems show a va- 
riety of physical conduction mechanisms. Due to their small size quantum mechanics becomes 
important and Coulomb blockade, interference, and Kondo physics are observed. Molecular 
systems fl]-^| are characterized by large electronic energies beyond room temperature and 
therefore offer the possibility to measure positions of molecular orbitals by transport spec- 
troscopy. Furthermore, mechanical degrees of freedom in molecules have energies of order 
1 — 10 meV which can be probed experimentally for temperatures in the Kelvin regime. This 
letter addresses the latter topic and combines nano-electronic with nano-mechanical proper- 
ties, in particular we model the experiment in Ref. Q . We find that local bosonic excitations 
have an important influence on single electron tunneling, and thus need to be included in 
models for molecular electronics. Based on those systems we propose another way to realize 
negative differential conductance (NDC). 

Attempts to model transport through molecular nano structures up to now focused more 
on the electronic structure pxj] . In this letter we go one step further and include molecular 
vibrations (or any other local bosonic excitation) as well. We remark that our approach is 
fundamentally different from models where the vibration serves as a shuttle for the electrons 
p| , p^| . We address the case where the vibrational frequency is several orders of magnitude 
larger than the frequency associated with tunneling events. Moreover the physics is dominated 
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Fig. 1 - Theoretical model system: A local system consisting of a quantum dot and a harmonic 
oscillator is coupled to two electronic reservoirs and a boson bath. 

by charging effects, therefore the single charge tunneling picture Jl3| is more appropriate than 
the scattering picture |l4j . 

Experiment. In their experiment Park and coworkers measured the current through a 
single Ceo molecule. Using a break junction technique they placed the molecule between two 
gold electrodes, to which it is weakly bound, i.e. the molecule sits in a Lennard- Jones-like 
potential on the gold surface. They could vary the source-drain and a gate voltage, that 
shifted the lowest unoccupied molecular orbital (LUMO). Thus they have built a Cgo single 
electron transistor (SET). We briefly summarize their findings. The I-V curve displayed a 
large conductance gap, which could be varied with the gate voltage. Once the current was 
flowing a series of equally spaced (u>q 5 meV) little jumps followed, which was attributed 
to the vibrational excitation of the whole molecule (center of mass motion) in the Lennard- 
Jones potential discussed above. More excitations at higher energies could also be found 
and could be attributed to e.g. intra-molecular vibrations. The experiment was conducted 
at T = 1.5 K = 0.13 meV, the charging energy exceeded 100 meV, and the separation of 
the LUMOs is assumed to be of similar magnitude. From the current one can deduce the 
tunneling broadening T to be of 0(0.1 /«eV). 

Theoretical Model. Our model describes the dynamics of a local system, consisting of the 
electronic states and a single-mode bosonic excitation due to the harmonic potential, which 
accounts for the vibrations. The local system is coupled to two electronic reservoirs, a boson 
bath, which allows for dissipation and a gate that can shift the position of the MOs, see Fig. [j]. 
Hence the Hamiltonian reads 

H = H eTes + H p h + H mo i + i/ m ol-ph + #mol-eres- (1) 

The electronic reservoirs are described by non-interacting electrons with a constant density 
of states po, hence H crcs = ^ feoT e kara\. crr akar, where r € {L,R} labels the left and right 
electrode and a the spin. The bosonic reservoir is given by i? p h = ^2 q u q d q d q . The charging 
energy and the separations of the MOs is large enough to take the LUMO as the only par- 
ticipating MO. For the following discussion we take the LUMO to be spin-degenerate, and 
include the Coulomb interaction in a Hubbard-like charging term. With a projection operator 
for the vibrational states Pw — \l) we have 

#moi = € °( V a) c l c ° + Un m + J2 4 lh Pu- (2) 
The gate electrode is hidden in the dependence of e CT on the gate's potential. The energy levels 
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of the oscillator are ej lh — (I + 1/2)ujq. The coupling to the boson bath ensures relaxation and 
dissipation of the vibrational energy in a natural way, and is generically written as 

#mol-ph = gWqPll' (dq + d q ) . (3) 

qW 

For simplicity, the constants gw q are assumed to be independent of Z, I' and we introduce the 
coupling a via the spectral density a\uj\ s — 27r^ 9 |<7 9 | 2 <5(ti; — e q ). The exponent s depends on 
the underlying microscopic realization of the boson bath, i.e. acoustical phonons, fluctuations 
in the electro-magnetic environment, etc. In the spirit of the Caldeira-Legett model |la|, we 
assume from now on an Ohmic bath, i.e. s = 1. Finally the tunneling term, which includes 
the possibility to excite the whole Ceo molecule to a vibration against the surface during a 
tunneling event, 

tfmol-crcs = ^ ^ T W P H' ( a Lr C ° + h ' C ) • ( 4 ) 
harlV 

Here, T r describes the tunnel matrix element in the absence of vibrational states, and the 
factor tui accounts for vibrational excitations. Their microscopic origin can be understood Q 
similar to what is known as the Franck-Condon principle: The electron tunnels faster than 
the molecule can relax into a new state. When one electron tunnels onto the molecule, it is 
likely that the charge distribution around it, i.e. the image charges in the gates which serve 
as capacitors, is not uniform. Consequently there will be a force acting upon the charged 
Ceo- Hence the potential minimum xq is shifted slightly by a distance S. Therefore we 
need the matrix elements between the states tp Xo ,i' in the potential centered at Xo and the 
states ip Xo +6,l of the potential shifted by the distance 5, i.e. tu> ~ (ipx +S,l\ipx ,l') ■ F° r the 
system under consideration the states are well described by the harmonic oscillator ones. 
Obviously it gives 6i t i' for S = or zero for S large, and we find that only a few transitions 
have considerable weight. With the tunneling broadening r r = 2irpo\T r \ 2 , we introduce the 
couplings r[ z , = r r |rn' 2 and Tu> — J2 r Y\ v , which control both the current and the excitations. 

Theory. The dynamics of the LUMO and all local bosonic excitations can be captured 
by a kinetic equation for the reduced density matrix, which is obtained by tracing out the 
reservoir and bath degrees of freedom of the full density matrix. For the parameters of interest 
double occupancy and off-diagonal matrix elements can be neglected (Eq ^> {T, r,w } and 
o^o S> r). Then only the rates describing electron tunneling including vibrational excitation, 
and those related to bosonic transitions are required. Note that they are completely distinct, 
since the charge on the molecule changes by one for the electronic ones, and remains constant 
for the boson-mediated ones. For weak tunneling r -C T and weak dissipation a <C 1, we can 
perform perturbation theory in the coupling constant in both cases. The rate describing a 
transition from zero charge and vibrational state V to charge one with spin a and vibrational 
state I is given by £ CT j,oj' = J2 r ^ai,M' witn 

Ki^^^fr^ + e^-e^). (5) 

For the reverse process the Fermi function f r (uS) will be replaced by (1 — f r (ui)). The rate 
describing phonon absorption and an accompanying vibrational excitation from I' to I (without 
change of the electronic state) is given by 

X+,=a\l'-l\u N(e] ib -eJ, ib ). (6) 

For the emission rate the Bose function N(u>) is replaced by 1 -I- N(ui). The kinetic equation 
for the probabilities of the state s reads P s = = (S ss ' P s ' — S S / S P S ), and the current is 
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obtained via 

1 = e E ( P « S ^',o ; - ^I'EqU) . (7) 

Results and Discussion. Before we focus on the specifics we discuss the general properties 
of the I-V characteristics. We tune the gate voltage such that the LUMO is above the Fermi 
levels. The voltage is applied symmetrically, i.e. fiL.R = ±eV/2. For low boson temperatures 
(solid curve in Fig. |J), absorption of phonons is suppressed and spontaneous emission, which 
in many molecules happens on a shorter time scale than the electron tunneling []l6f , will force 
the molecule to stay predominantly in its vibrational ground state. In this case, as in the 
experiment, we see no current until eVbias — followed by a series of small steps with 
distance 2oj . The latter are due to excitations of higher vibrational states while tunneling. 
Each time the bias voltage exceeds 2(e CT + nuj ), a new tunneling-in process |0, 0) — ► \a, n) can 
occur leading to an increase of the current. The excited vibrational state relaxes immediately 
to the ground state n = and tunneling-off happens via the process \a, 0} — > |0,0). The 
number of steps is determined by the number of vibrational states taken into account. The 
width of the rise for each single step is controlled by max(T c i,r), which in our case is T e \. In 
Fig. H we show results where all vibrational states can be accessed equally, i.e. tu> = 1. We 
remark that our results for strong dissipation are not special for the Ohmic bath, but rather 
independent of (positive) s, since the energy is dissipated immediately. 

The situation is more complicated in the absence of dissipation (dashed curve in Fig. ||). 
Here, once the threshold for the onset of current has been reached, the tunneling-off processes 
|cr, 0} — > |0, n), with ncoo < e CT , can leave the initially not excited molecule in an excitated state. 
This state then serves as a starting point for other tunneling-in processes 0, n) — > \a, to), with 
to < n, leading to an increase of the current compared to the case with dissipation (after each 
additional step in the I-V-characteristics, the allowed values for to increase by one). 

For high boson temperatures where emission and absorption of phonons is possible (dot- 
dashed curve in Fig. |), the I-V curve for eVbias > is similiar to the case without dissipation. 
However, below the treshold, current precursors occur due to tunneling-in processes where the 
vibrational state lowers its energy. The higher vibrational states have here a finite probability 
due to absorption of phonons which is not possible in the absence of dissipation. However, 
since this is not seen in the experiment, one can rule out such high boson temperatures. 

Experimentally the number of steps that can be resolved is much less than 10 and we there- 
fore proceed by taking into account the dependence of the Franck-Condon factor tu> on the 
vibrational states. Assuming that the potential is well described in a harmonic approximation 
with parameters from the Lennard-Jones potential jTffl, then only the shift <5 is unknown. A 
crude estimate, based on the force due to an image charge in the electrode, is given by Park 
et al^. Going beyond this picture is difficult because the electrodes' geometry is not known. 
However, our results show that a better estimate may be obtained. In Fig. || we show the I-V 
characteristics for three 5's of the order of the previous estimate of 4 pm. The ratio of the first 
and second step height is highly sensitive to S. Comparing with the experiment one therefore 
would suggest a value of S w 3 pm. Moreover it can be seen that the current saturates much 
faster than without Franck-Condon factors. This is not surprising, because the molecule is 
mostly in its ground state from which the higher states can not be reached. This is different 
without dissipation, where even a small NDC can be realized. The differential conductance, 
which shows a peak for every excitation, is shown in Fig. [|, with the NDC effect in the inset. 
We remark that our results agree very well with the experimental data ||. However, there 
remain features which can not be reproduced in our calculation, namely a rising background 
and more excitations at higher energies. The origin of the background is not precisely known 
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Fig. 2 - Influence of the boson bath on the I-V characteristics: Emission only (solid, T p h = 0.13 meV), 
no phonon bath (dashed), and emission and absorption (dot-dashed, T p h = 13 meV). T e l — 0.13 meV, 
e CT = 15 meV, a = 0.05, uo = 5 meV, r r = 0.015 fieV, 10 vibrational states. Transitions are indicated 
for the solid and dashed line. 



yet, the other excitations will be addressed below. 

The NDC effect as seen before is due to a new mechanism specific for additional bosonic 
degrees of freedom. To exhibit the physics, Fig. ^ shows the effect for a simplified model of one 
electronic level with two vibrational states I = 0,1. The T = currents of the two plateaus 



0.2 r 



0.15 115 



% 0.1- 



0.05 





Fig. 3 




\r ■ 

I 
I 



20 



40 



60 

[meV] 



80 



CM 

CD 



> 





I ' le-05 








5e-06 


\ — n — <— 















50 60 70 

s 


J 


! i 



100 



30 



Fig. 4 



40 



bias 



50 

[meV] 



Oil 



70 



Fig. 3 - Including the Franck-Condon factors with 5 = 3 pm (solid), 8 = 4 pm (dashed), and 8 — 5 pm 
(dot-dashed). Inset: Comparison of with (lower) and without Franck-Condon factors (upper curve). 
T c i = Tph — 0.13meV, e CT = 15 meV, a = 0.05, u>o = 5 meV, T r = 0.1 neV, 10 vibrational states. 

Fig. 4 - dl-dV curve, 8 = 3 pm with bath (solid) and without (dashed in inset), and 5 = 4 pm with 



bath (dot-dashed). T e i = T ph = 0.13 meV e a = 15 meV, a = 0.05, lu = 5 meV, F u 
vibrational states 
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Fig. 5 - Negative differential conductance with two vibrational states. Main panel shows Tn ranging 
from 0. Of Too (bottom curve) viaO.f Too, 0.2r 00 , 0.5r O o, 0.7r 00 , 0.8r o0 to 0.9r o (top curve), a = 0. 
fnset: a = 5 f 0~ 5 (top) down 5f0 -9 . Each line corresponds to one order of magnitude. Common 
parameters are T e \ = T ph = 0.f3 meV, e a = f5 meV, loo — 5 meV, Toi = Tio = O.f Too and T r = O.f 
^eV. 

Fig. 6 - Negative differential conductance with TV vibrational states and r»/ = — l'\ + f ) 2 f /((I + 
l')/2) 2 . Parameters are T ei = T ph = O.f 3 meV, e CT = 2 meV, lu = 5 meV, and T r = O.f ^eV. 



are (for T i = T 



10) 



v r i t) /"p \ o-n r 2r oi(roo + r i + 2rn) + rn(3r o + Tn' 

-TOO 1 00 + "0l( i 01 + 1 11 J — -fill 01 — 1 oo — 



6r 01 roo + 2r 11 (r i + 4roo) 



hn - Poo(r 00 + r i) + P i(r i + r u ) = x - (r 00 + 2r 01 + r n ) . (8) 

In the limit Toi Too we see that only Tu < Too is needed to observe NDC. For Too 3> 
Tio = Toi = Tu, we get /fi n /^mid = 14/20. The origin is that for the second plateau the 
probability distribution changes such that even weakly coupled vibrational states are occupied 
equivalently to the strongly coupled ones. However, their contribution to the current is small, 
and therefore the current decreases. We stress that the Tu> naturally show the appropriate 
kind of behavior. Assuming a simple functional form for them we show in Fig. ^ how many 
bosonic levels can in principle lead to large peak to valley ratios. For multi-mode bosonic 
excitations the effect survives and can even be more drastic, when many excitations are energy 
degenerate. NDC has already been seen in other nanostructures flg|-|2"o(|, and explanations in 
terms of spin-blockade |^lj and asymmetric coupling |^2| have been suggested. In our case 
however, no additional electronic states at higher energies are required, nor do we need the 
full many-body state of the structure. 

Finally we comment on the influence of bosonic excitations in general molecular electronic 
devices. Multi-mode bosonic excitations can be included in the same way. Other degrees of 
freedom could be due to internal vibrations, e.g. the torsional vibration of the oligo-phenyl 
molecules used in the experiment by Chen and coworkers |2C| ] have energies in the meV range 
and should be observable at low enough temperatures. The structure of smaller bosonic 
excitations would display itself on top of the larger ones, provided they are well separated. 
For heavier molecules where uj ~ T we expect for T < T an interesting interplay including 
interference effects between the various excitations, partly related to off-diagonal elements 
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of the density matrix. In addition, for T ^> T, this can lead to an increasing background 
vibrational noise. This background (not in the gap) is indeed seen in the experiments and 
it is quite likely that such small local bosonic excitations play an essential role in it. We 
remark that this kind of behavior is fundamentally different from artificial semiconductor 
nano structures, which are embedded in a substrate. 

Summary. We described single electron tunneling through a system with local bosonic 
degrees of freedom and dissipation to external heat baths. We applied our model to a Can 
molecule, which can be excited to center of mass oscillations. We found that our results 
compare well with the experiment. We showed that this model may be used to realize negative 
differential conductance. We concluded that transport through molecular nano-devices is 
inevitably controlled by its electronic and mechanical structure. 
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